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New  results  are  presented  for  data  traffic  performance  in  an 
integrated  multiplex  structure  which  Includes  circuit-switching  for 
voice  and  packet-switching  for  data.  The  results  are  developed  both 
through  simulation  and  analysis,  and  show  that  excessive  data  queues  and 
delays  will  build  up  under  heavy  loading  conditions.  These  large  data 
delays  occur  during  periods  of  time  when  the  voice  traffic  load  through 
the  multiplexer  exceeds  its  statistical  average.  A variety  of  flow 
control  mechanisms  to  reduce  data  packet  delays  are  investigated.  These 
mechanisms  include  control  of  voice  bit  rate,  limitation  of  the  data 
buffer,  and  combinations  of  voice  rate  and  data  buffer  control.  Simu- 
lations indicate  that  these  flow  control  mechanisms  are  quite  effective 
in  improving  system  performance.  A combination  of  data  buffer  limitation 
with  data-queue-dependent  voice  rate  control  was  the  most  effective  flow 
control  technique  tested. 
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INTRODUCTION 


I . 


An  analysis  of  the  performance  of  an  Intagrated,  circuit-  and 

packet -switched  multiplexer  structure  was  presented  recently  by  Fischer 

l 2 

and  Harris.  The  structure,  Introduced  by  Covlello  and  Vena,  is  based 

on  a time-slotted  frame  format  where  a certain  portion  of  the  frame  is 
allocated  to  ctrcutt-switched  calls  and  the  remaining  capacity  is 
reserved  tor  data  packets.  To  increase  channel  utilisation,  a "movable 
boundary"  feature  is  included,  so  that  data  packets  are  allowed  to  use 
any  residual  c trcuit-switched  capacity  momentarily  available  due  to 
voice  traffic  variations.  The  performance  measures  employed  in  the 
analysis  wore  blocking  probability  for  circuit-switched  voice  customers 
and  average  waiting  time  for  data  packets.  Results  indicated  signi- 
ficant saving,  in  terms  of  channel  requirements  for  the  integrated 
system,  as  compared  to  a system  with  separate  circuit-  and  packet- 
switched  facilities  handling  the  same  traffic  loads  with  the  same 
blocking  probability  and  average  packet  waiting  time.  However,  an  error 
has  been  discovered  in  the  original  derivation  of  average  data  packet 
waiting  time.  New  results'  indicate  much  larger  data  packet  delays  than 
previously  predicted  for  the  case  where  the  service  capacity  reserved 
for  data  is  smaller  than  the  data  traffic  load.  These  results  tend  to 
indicate  that  the  advantage  of  the  integrated  system  in  terms  of  channel 
utilization  Is  much  loss  than  indicated  previously. 

However,  these  results  should  certainly  not  he  considered  in 
isolation  in  evaluating  Integrated  system  approaches.  In  particular, 
the  analysis  here  assumes  constant  levels  of  average  voice  and  data 
traffic  load,  and  shows  that  a dynamically  movable  boundary  is  not  very 
effective  in  allowing  data  packets  to  utilize  free  capacity  due  to 
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momentary  fluctuations  in  voice  traffic.  However,  it  is  also  shown  In 
this  paper  that  appropriate  t low  control  mechanisms  can  appreciably 
Improve  data  traffic  performance  with  a dynamically  movable  bounJary. 

In  addition,  an  integrated  system  ought  to  allow  effective  adaptation  to 
the  more  gradual  variations  in  voice  and  data  traffic  loads  that  would 
typically  occur  over  a 24-hour  period,  finally,  other  criteria  such  as 
flexibility  and  overall  hardware  cost  must  be  considered  together  with 
efficiency  of  channel  utilization  in  comparing  various  system  configurations. 

The  primary  purpose  of  this  paper  Is  to  develop  modified  results 
for  the  data  performance  of  the  multiplexer  structure.  These  new 
results  have  been  obtained  through  two  approaches:  (1)  a simulation 
designed  to  match  the  original  model  as  closely  as  possible;  (2)  an 
analysis  of  a somewhat  modified  system  model,  which  Is  more  analytically 
tractable  yet  retains  the  same  essential  performance  characteristics. 

The  simulation  model  and  sample  results  will  be  presented  first,  and  an 
attempt  will  be  made  to  explain  the  basic  cause  of  the  largo  packet 
delays.  This  material  is  contained  in  Section  II.  Then  the  analytical 
approach  together  with  some  sample  results  Is  given  in  Section  III.  The 
second  purpose  of  the  paper  Is  the  Investigation  of  techniques  for 
controlling  the  flow  of  voice  and/or  data  traffic  through  the  multiplex 
structure,  with  the  goal  of  decreasing  data  packet  delay.  This  work  is 
described  in  Section  IV.  Explanation  of  the  error  in  the  earlier 
analysis  and  a more  detailed  description  of  the  new  analytical  approach 
as  well  as  comparisons  with  the  simulation  runs  are  given  in  Appendices. 


II.  SIMULATION  OF  MULTIPLEXER  PERFORMANCE 


A . Frame  Structure 

The  multiplexing  technique  under  Investigation  is  a Time  Division 
Multiplex  (TDM)  scheme  in  which  time  slices  of  fixed  duration,  a frame 
period,  are  partitioned  and  allocated  to  the  transmission  of  digitized 
, voice  and  data  packets.  The  voice  component  of  a frame  is  further 

divided  into  slots  allocated  to  ongoing  line  switching  communications. 

' In  general  the  voice  slots  can  he  of  varying  size  depending  on  the  bit 

rates  of  different  voice  users,  and  the  size  of  data  packets  is  also 
variable.  For  simplicity,  and  in  order  to  allow  direct  comparison  with 
the  earlier  analysis, ^ the  simulation  here  has  assumed  a multiplexer 
frame  (taken  to  be  of  duration  b - 10  ms)  divided  into  S + N equal- 
capacitv  slots.  The  nominal  bit  rate  for  the  voice  coders  is  taken  as 
8 kbps,  which  is  accommodated  with  an  80  bit  slot  size  and  a 10  ms 
frame  period.  Voice  traffic  is  allowed  to  occupy  the  first  S slots,  and 
the  remaining  N slots  are  reserved  for  data  packets  which  arc  each 
assumed  to  occupy  one  slot.  A voice  call  retains  its  slot  for  the 
duration  of  the  connection;  but  if  all  S voice  slots  are  busy  when  a 
new  call  is  initiated,  that  call  is  blocked.  Data  traffic  is  permitted 
to  use  any  voice  slots  which  may  be  momentarily  free  due  to  statistical 
fluctuations  in  the  voice  traffic. 

B.  Voice  Traffic  Simulation  Model 

Voice  traffic  is  modelled  by  a Poisson  call  arrival  process  of 
rate  \ and  exponentially-distributed  call  holding  times  with  mean  p *• 

The  average  number  of  call  arrivals,  \b,  in  a frame  interval  is  generally 
small  enough  so  that  the  frame  structure  can  be  ignored  for  voice 
traffic  analysis.  Thus  the  voice  multiplexer  can  be  described  as  a 
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classical  Harkov  S-server  loss  system  M/M/S/S  whoa*  state-transit  Ion- 
rate  diagram  is  ahown  in  Fig.  1. 

Let  ttyCt)  represent  the  number  of  calls  in  progress  at  a given 
time.  A simulation  of  the  variation  of  n (t)  with  time,  based  on  the 
state-transition  diagram,  has  been  developed  as  follows.  Assume  that 
nv(t)"k  at  a particular  starting  time.  Then  nylt)  Is  held  at  k for  a 
time  t drawn  from  an  exponential  distribution 

fx  (t)  - Tk"1  e'l/lk,  (1) 

where  the  mean  holding  time  Is  determined  as 

Tk  - 1/ (\  + kn)  k - 0,  1,....  S-l  (2) 

T j,  - 1/Svi. 


After  a time  f,  «y(t)  is  increased  to  k-H  with  probability 


P (k)  - A/ (X  + ku)  k - 0,  1 S-l  O) 

up 

P (S)  - 0 
up 


or  decreased  to  k-l  with  probability  1-P  fk).  This  process  Is  repeated 
as  often  as  desired  to  generate  sample  functions  of  nyftX. 

It  is  useful  to  note  that  the  steady  state  blocking  probability  for 
the  S-server  loss  system  Is  given  by  the  "Krlang-h  formula" 


pt . 

i (X/L)*Vk, 

k-0 


f4a> 


4 


and  that  the  expected  value  of  nv><nv>,  Is  given  a* 


<ny>  - (X/p)  (1-PL), 


(4b) 


where  ny  la  the  random  variable  representing  the  number  of 
voice  calls  In  the  system. 

C.  Data  Traffic  Simulation  Model 

Data  packets  are  assumed  ;o  arrive  in  a Poisson  process  with  rate  0 

packets/ sec.  Packet  size  Is  assumed  fixed  and  equal  to  the  number  of 

bits  In  one  slot.  During  any  frame,  the  number  of  slots  available  for 

transmission  of  data  packets  Is  N+S-n^(t),  the  sum  of  N dedicated  data 

slots  and  S-n  (t)  unused  voice  slots, 
v 

Consider  a period  of  time  of  duration  t^  between  the  1th  and  (1+1 )st 
transition  In  nv(t),  during  which  n^it)  “ nv^;  the  number,  n+,  of  data 
packets  arriving  during  this  interval  Is  drawn  from  a Poisson  distribution 


P(n+) 


(0t1)n+  e~0ti 


‘V 


(5) 


with  mean  Ot^.  The  maximum  number  of  data  packets  which  can  be  sent  out 
during  this  time  is 


n_  - (N  + S - ny(1))  (t ±/b) 


(6) 


(since  generally  t^  » b,  the  fact  that  t^/b  is  not  necessarily  an 
integer  has  been  ignored.)  Let  the  number  of  data  packets  waiting  for 
service  at  the  beginning  of  the  ith  interval  be  n^^\  The  evolution  of 
this  random  variable  to  the  next  interval  is  simulated  as 


- max  (n^  + n - n , 0) 
a a + 


where  the  max  insures  that  n^  ^ never  becomes  negative.  This  equation 


is  strictly  valid  as  long  as  the  data  queue  does  not  empty  during  the 


interval  t^.  If  the  data  queue  does  empty  during  the  interval,  then 


may  be  slightly  greater  than  predicted  by  (7),  if  the  dispersion 


of  packet  arrival  times  causes  some  packet  service  slots  to  go  unused. 


This  slight  discrepancy  during  periods  when  n,  is  close  to  zero  would 

d 


have  negligible  effects  on  the  results  presented  below,  and  has  been 


ignored. 


The  model  used  for  the  combined  voice/data  simulation  is  described 


by  equations  (l)-(7).  Note  that  nd  is  updated  only  at  times  of  transi- 


tion in  n (t).  For  the  purpose  of  display  and  of  determining  time 


averages,  nj(f)  is  assumed  to  vary  linearly  between  its  sample  values 


n^  . Suitable  random  number  generators  are  used  to  determine  T^,  n+. 


and  the  up/down  transition  decisions  for  nv(t).  The  average  data  packet 


waiting  time  is  determined  from  measured  values  of  <nd>  by  Little's 


theorem  as 


Wd  = <nd>/0  , 


where  <nd>  is  a time  average  of  nd> 


D.  Simulation  Results 


1.  Analysis  of  Sample  Runs 


Sample  functions  of  nv(t)  and  nd(t)  are  shown  in  Figure  2 for 


S = 10,  N = 5.  The  average  call  holding  time  is  taken  as  1/y  = 100  sec. 


with  X = .05  sec  so  that  the  offered  voice  traffic  is  X/y  Erlangs. 
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TIME  (sec) 


Fig.  2.  Samole  functions  of  n (t)  and  n,(t)  for  S = 10,  N = 5, 

X « 0.05  sec  , y ■ 0.01  sec  1 and  0 = 900  packets/sec  (0b  = 9 packets/ 


Application  of  (Aa)  yields  • .018  and  <nv>  ■ A. 91  so  that  the  average 
number  of  alota  available  for  data  packets  la 

n~  -N+S-<n>-  10.09. 
p v 

The  data  packet  arrival  rate  la  0 - 900  aec  * for  an  average  data  utili- 
zation factor  of 

- 0b/np  - .89. 

The  plots  represent  2500  sec  of  real  time,  during  which  approximately 
230  voice  calls  either  entered  or  left  the  system.  It  Is  clear  that 
nv(t)  is  a highly  correlated  (in  time)  random  process  which  exhibits 
swings  of  long  duration  above  and  below  <nv>.  The  error  (see  appendix 
Al)  in  the  previous  analysis  effectively  caused  the  time  correlation  of 
n^(t)  to  be  neglected,  leading  to  large  underestimates  of  data  packet 
waiting  time. 

During  "idle"  periods  where  nv(t)  is  low  enough  so  that  0b  < N + 

S - nv(t)  more  than  enough  capacity  is  available  to  handle  incoming  data 
traffic  and  an  initially  empty  data  queue  nd(t)  will  remain  essentially 
empty.  But  will  build  up  significantly  during  busy  periods  where 

Ob  > N + S-nv(t).  For  this  example,  a long  busy  period  during  which 
nv(t)  > 6 for  about  500  sec  begins  at  about  t “ 900  sec.  During  this 
period  the  data  queue  builds  up  to  about  30,000  packets.  The  data  queue 
is  eventually  emptied  during  a subsequent  idle  period  of  the  voice 
channel,  but  then  builds  up  again  during  the  next  busy  period.  The 
average  buffer  size  during  a run  of  which  the  first  half  is  depicted  in 
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r — • 1 1 

Fig.  2 is  ■ 9200  packets,  corresponding  to  a mean  delay  of  - 

10.2  sec. 

The  behavior  depicted  in  Fig.  2 is  typical  of  that  exhibited  in 

many  similar  runs.  A segment  of  another  pair  of  sample  functions  nv(t), 

nj(t)  is  shown  in  Fig.  3 for  S ■ 50,  N ■ 25,  X - .4,  U ■ .01,  0 ■ 3400. 

Here  <n  > - 39.2  and  p,  - .02.  There  are  S-<n  > - 10.8  residual  voice 
v L v 

channels  on  the  average  in  steady  atata,  and  the  average  data  load  fills 
9 of  these  residual  channels.  The  plots  cover  500  sec  of  real  time  and 
about  390  voice  transitions.  The  long  busy  period  ending  near  t ■ 400 
leads  to  a maximum  data  buffer  of  about  63,000  packets.  Average  buffer 
sire  for  this  run  is  <nj>  ■ 12,852  packets,  and  average  delay  is  - 
<n^>/8  s 3.78  sec. 

Figure  4 illustrates  what  happens  in  the  extreme  case  when  only  one 
slot  per  frame  (S-l,  N-0)  is  available.  Here  X - .01,  p ■ .01,  and 
<nv>  - .5  so  that  the  slot  is  occupied  by  voice  half  the  time.  The  data 
load  8b  is  .4  packets/slot  so  that 

p,  - 8b/n  ■ .4/. 5 " .8. 
d p 

The  data  buffer  n.(t)  obviously  must  increase  while  n (t)  “ 1 and  these 
d v 

increases  must  be  worked  off  while  nv(t)  “ 0.  For  this  example, 

<n^>  " 7047  packets  and  • 176  sec.  The  effect  of  the  average  voice 
holding  time  1/p  should  be  apparent  from  this  example.  If  1/u  were 
doubled  to  200  sec  with  X/p  fixed,  a typical  sample  function  of  nv(t) 
would  be  as  shown  in  Fig.  4,  except  that  twice  as  much  elapsed  time 
would  be  represented.  With  8 fixed,  this  time  scale  change  would  lead 
to  approximately  a doubling  of  the  values  plotted  in  Fig.  4,  and 
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- 40,000 

M 

l 

' 20,000 


Fig.  4.  Samnle  i 
X » 0.01  sec  , U 
f rame) . 


therefore  a doubling  of  <n^>  and  W^.  The  Increased  Is  an  obvious 
consequence  of  the  fact  that  data  packets  (arriving  when  nv(t)*l)  must 
wait  twice  as  long  on  the  average  for  the  ongoing  conversation  to 
terminate. 

Changes  in  p for  fixed  A/p  will  have  similar  effects  on  the  results 
depicted  in  Figures  2 and  3.  Actually  the  holding  time  1/p  * 100  sec  is 
shorter  than  typically  encountered  in  telephone  traffic,  so  that  the 
results  given  here  on  data  delays  are  probably  optimistic  by  about  a 
factor  of  2 or  3. 

2.  Average  Data  Performance  Statistics 

Average  data  packet  delays  as  measured  by  the  simulation  are 

depicted  In  Fig.  5 for  S“10,  N«5,  \/p  « 5,  and  p ■ .01.  Each  plotted 

point  represents  an  average  of  4 runs,  with  each  run  comprising  5000 

transitions  in  nv(t),  or  about  50,000  sac  of  real  time.  Actually,  the 

average  buffer  size  ^nj>  was  measured,  and  (8)  used  to  obtain  W^.  For 

example,  at  0b  “ 9,  1 10  sec  and  <n^>  **  9000  packets.  In  each  run, 

the  maximum  as  well  as  the  average  of  n^(t)  was  measured.  Cenerally 

(n.(t))  » <n,>,  as  indicated  in  Table  1.  (Packet  size  has  been 

d max  d 

taken  as  80  bits  as  in  the  previous  analysis.*)  Enough  storage  must  be 

allocated  to  hrndle  (n.(t))  . For  example,  with  Ob  - 9.0,  about 

d max  r 

800,000  words  of  storage  would  be  needed  to  prevent  data  buffer  overflow. 

Results  were  presented  previously*  for  data  waiting  times  for  the 
same  values  of  N,  S,  and  \/p  as  utilised  in  obtaining  Fig.  5.  Assuming 
b - 10  ms,  the  predicted  waiting  time  was  ® 20  ms  for  Ob  - 9.  This 
predicted  delay  is  about  a factor  of  500  below  the  simulation  results 
obtained  here. 
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5 6 7 

fib  (pockets /from* ) 

Fig.  5.  Average  data  packet  delays  as  function  of  data  load  for  S - 10, 
N - 5,  X ■ 0.05  sec  p ■ 0.01  sec-1.  Each  plotted  point  represents  an 
average  of  4 runs,  with  each  run  comprising  5000  transitions  in  n (t)  and 
about  50,000  sec  of  real  time. 
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Average  data  packet  delays  for  S - 50,  N - 25,  X/p  - 40.0  and 
■ .01  are  shown  in  Figure  6.  Delays  are  slightly  less  than  In  Fig.  5 
for  the  same  values  of  • 0b/ (N  + S - <n^>),  partially  because  the 
larger  value  of  A/p  Implies  that  transitions  In  nv(t)  occur  at  a faster 
rate.  However,  the  larger  values  of  9 Imply  that  buffer  sizes  are 
similar  to  those  shown  In  Table  1.  Again,  for  > .9,  average  waiting 
times  on  the  order  of  10  sec  are  exhibited. 


( sec) 


r 


Kig.t).  Average  data  packet  delays  as  function  of  data  load  for  S - ^0, 

N - 25,  X ■ 0.4  sec-1,  and  p " 0.01  sec-1.  Each  plotted  point  represents 

an  average  of  2 to  4 runs,  with  each  run  comprising  5000  transitions  in 

n (t>,  and  about  6500  sec  of  real  time, 
v 
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TABLE  1 

AVI: RACE  AND  MAXIMUM  PATA  BUFFER  SIZE. 

IN  THOUSANDS  OF  lb- BIT  WORDS , 

AS  FUNCTION  OF  PATA  LOAD 

(Each  entry  represents  an  average  of  sn,>  or  (n.l 

d d max 

over  4 runs,  with  each  run  comprising  5000  transitions 
in  n v ( d ) and  about  50,000  sec  of  real  time.) 


,,  1 

Ob 

(packets/ 
f rame I 

Average 

Data  Buffer 
(K  words) 

Max imum 
Data  Buffer 
(K  words) 

III.  MODIFIED  ANALYTICAL  MODEL 


A.  Model  and  Notation 

If  one  Ignores  the  effect  of  the  time  quantization  introduced  by 
the  frame  structure  and  assumes  exponentially-distributed  rather  than 
constant  data  packet  lengths,  the  entire  system  can  be  analyzed  as  a 
two-dimensional  Markov  chain.  It  is  useful  to  introduce  new  notation 
where  the  voice  customers  are  denoted  by  class  1 and  the  data  customers 
are  class  2.  Assume  the  arrival  process  for  each  is  Poisson  with  rate 
X^i  “ 1,2)  and  the  service  time  is  exponential  with  rate  p^.  Let  Qj 
denote  the  number  of  class  i customers  in  the  system  and  P represent 
the  steady-state  probability  that  - i and  Qj  ■ j.  Referring  to  the 
simulation  model  described  above,  X^  corresponds  to  X,  to  p,  X^  to 
0,  p_  to  b X,  Q,  to  n . and  Q~  to  n.. 

I IV  Z d 

There  are  S + N server  channels  with  voice  calls  having  priority  on 
S of  the  channels.  If  a voice  call  shows  up  and  < S,  the  call  will 
seize  a free  channel  if  there  is  one,  otherwise  it  will  preempt  a data 
packet  which  returns  to  the  head  of  the  data  queue.  No  queueing  is 
allowed  for  voice.  In  Appendix  II,  a general  set  of  two-dimensional 
difference  equations  for  P^  is  written,  and  a solution  procedure  is 
outlined.  The  performance  parameter  of  interest  is  the  expected  waiting 
time  for  data  packets,  which  is  related  to  P according  to 


W.  - T~  E[Q  ] - f-  E E j P . 
d X2  2 X2  i-0  J-0  13 


(9) 


B.  SinRle-Channel  Example 

For  the  special  case  of  one  voice  channel  and  no  data  channels 
(S  - 1,  N - 0),  a reasonably  simple  expression  for  emerges,  in  the  form 
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(10) 


P2d  + pj)2  + Px  ^2/p1 
Wd  “ X2(l  + p1)(l  - P2  - PjP2) 

where  p^  ■ X^/pj  anJ  p2  " ^2^2’  The  *actor  ^2^1  in  c^e  numerator 
represents  the  expected  number  of  data  packets  to  arrive  during  an 
average  voice  holding  time.  This  factor  can  be  quite  large  and  accounts 
in  large  part  for  the  very  large  data  delays  noted  in  the  simulation. 
Applying  (10)  to  the  case  illustrated  in  Fig.  4 where  S - 1,  N - 0, 

X^  ■■  .01,  pj  • .01,  X2  - 40,  and  u2  • 100  results  in  * 250  sec.  This 
analytic  result,  based  on  exponentially-distributed  packet  lengths,  is 
reasonably  consistent  with  the  result  of  the  simulation,  where  with 
constant  data  packet  lengths  was  estimated  at  176  sec. 
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IV.  INVESTIGATION  OF  FLOW  CONTROL  TECHNIQUES 


The  conclusion  that  can  be  drawn  from  the  simulation  and  analytic 
results  shown  above  Is  that  attempts  to  achieve  high  channel  utilization 
can  lead  to  data  packet  queues  so  large  as  to  overflow  any  reasonable 
amount  of  storage  in  the  multiplexer.  Even  assuming  infinite  storage, 
the  mean  data  delays  are  large  because  of  the  buffer  buildup  during 
periods  when  voice  channel  occupancy  is  high.  The  need  for  a flow- 
control  mechanism  is  apparent,  and  two  types  of  flow  control  have  been 
investigated:  (1)  voice  rate  control,  where  voice  coders  are  assumed  to 

be  capable  of  operating  at  a variety  of  rates  and  a new  voice  call  is 
assigned  (at  dial-up)  one  of  several  available  bit  rates,  based  on  the 
current  voice  utilization  and/or  data  queue  size;  and  (2)  data  flow 
control,  where  a fixed  limit  is  imposed  on  the  size  of  the  data  queue. 

A set  of  simulations  has  been  run  to  gather  statistics  on  multi- 
plexer performance  with  different  flow-control  strategies.  Generally 
total  channel  capacity  was  fixed  at  120  kbps,  with  voice  capacity 
limited  to  80  (sufficient  to  accommodate  N « 10  users  at  the  nominal  8 
kbps  rate),  and  40  kbps  dedicated  to  data.  For  comparison,  one  of  the 
flow  control  schemes  was  tested  with  600  kbps  channel.  A movable 
boundary  was  employed  to  allow  data  packets  to  be  sent  in  temporarily 
unused  portions  of  the  voice  capacity.  The  voice  traffic  was  modelled 
as  above,  with  X set  at  .05  sec  ^ and  p ■ .01  sec  \ corresponding  to  5 
Erlangs  of  offered  voice  traffic.  The  maximum  number  of  simultaneous 
voice  calls  was  limited  to  10,  with  blocked  calls  cleared.  Data  packets, 
with  fixed  lengths  of  80  bits,  were  assumed  to  arrive  in  a Poisson 
process.  The  purpose  of  the  simulations  was  to  Investigate  system 
performance  for  data  packet  arrival  rates  greater  than  500  per  second 


(or  equivalently  5.0  per  frame),  when  a portion  of  the  voice  capacity 
must  be  utilized  for  data.  For  the  cases  where  the  data  packet  arrival 
rate  Is  less  than  500  per  second  no  large  data  queues  build  up.  The 
performance  measures  of  interest  for  data  include  average  delay,  average 
queue  size,  maximum  queue  size,  and  fraction  of  packets  arriving  when 
the  data  queue  is  full  (when  data  flow  control  is  employed).  With  voice 
rate  control,  the  distribution  of  callers  among  the  available  bit  rates 
and  the  average  voice  bit  rate  assigned  serve  as  performance  measures 
for  voice  users. 

A.  Experiments  with  Voice  Rate  Control  but  No  Data  Flow  Control 

Up  to  now,  all  voice  users  were  assumed  to  operate  at  8 kbps.  The 
voice  channel  utilization  exhibited  large  variations  around  its  mean, 
and  data  queues  would  build  up  during  the  peaks  of  voice  utilization. 

The  idea  behind  the  voice  rate  control  techniques  studies  here  was  to 
cut  down  the  peaks  of  voice  channel  utilization  by  assigning  lower  bit 
rates  to  callers  who  enter  the  system  when  utilization  is  high.  This 
scheme  of  course  assumes  that  each  voice  user  has  a flexible  vocoder 
capable  of  a variety  of  rates.  The  presumption  is  that  the  assigned 
voice  bit  rate  and  the  corresponding  speech  coder  performance  will  vary 
depending  on  the  traffic  load  in  the  system. 

The  results  for  four  voice  rate  control  schemes  are  reported  here. 
To  define  these  schemes,  the  following  notation  is  introduced: 

Ry  = sum  of  bit  rates  of  all  voice  users  active  at  a given  time; 

Cy  * maximum  channel  capacity  allocated  for  voice; 

" number  of  packets  currently  in  the  data  queue; 

R * bit  rate  assigned  to  a new  voice  caller  at  dial-up. 
new  r 

Unless  otherwise  noted,  Cy  was  fixed  at  80  kbps.  The  first  voice  rate 
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control  scheme  that  was  tried  is  illustrated  in  Fig.  7.  New  calls  were 
assigned  at  2,  4,  8 or  16  kbps  as  a function  of  current  total  voice 
channel  utilization  Ky/Cy,  with  lower  rates  assigned  during  higher 
utilization  periods.  The  inclusion  of  a 16  kbps  rate  will  provide 
better  voice  service  for  some  users,  but  has  the  effect  of  filling  in 
the  valleys  of  voice  utilization  that  are  present  when  voice  rate 
control  was  not  employed.  The  second  scheme  was  similar  to  the  first 
except  that  the  voice  rate  of  16  kbps  was  not  used,  but  8 kbps  was 

| 

assigned  to  callers  arriving  when  Ry/Cy  was  less  than  .25.  The  third 
approach  used  voice  rate  control  combined  with  monitoring  of  the  data 
queue  was  included.  When  the  data  queue  was  empty,  voice  rates  were 
assigned  as  in  Fig.  7;  however  voice  callers  arriving  when  the  data 
queue  was  not  empty  were  assigned  a rate  of  2 kbps.  In  the  fourth 
technique,  the  rate  of  a new  call  was  governed  solely  by  the  size  of  the 
data  queue.  If  the  queue  was  empty,  a rate  of  8 kbps  was  assigned.  If 
the  queue  was  non-empty  but  did  not  exceed  150  packets,  4 kbps  was 
assigned.  Otherwise,  2 kbps  was  used.  The  four  approaches  just  out- 
lined will  be  referred  to  below  as  schemes  vl,  v2,  v3  and  v4;  vO  will 
denote  the  situation  where  no  control  is  imposed.  These  voice-rate 
control  schemes  are  summarized  as  follows: 


vO: 

R 

new 

■ 

constant  “ 8 kbps 

16  kbps 

Vcv  - ’25 

vl: 

R 

8 kbps 

•25  < Rv/Cv  < .5 

new 

4 kbps 

.5  < Rv/Cy  < .75 

2 kbps 

.75  < Rv/Cv  < l.i 

8 kbps 

Vcv  i -5 

v2: 

R 

new 

4 kbps 

•5  < VCV  1 -75 

.2  kbps 

.75  < Ry/Cv  < l.i 

CURRENT  VOICE  CHANNEL  UTILIZATION  Ry/Cv 

Mg. 7.  Kxample  ol  voice-rate  control  baaed  on  current  voice  channel 
utilization.  The  hit  rate  assigned  to  a new  user  is  plotted  as  a 
function  of  the  voice  channel  utilization  Ry/Cy. 
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v3:  R 


v4:  R " < 4 kbps 
new 

2 kbps 


| Assigned  as  in  vl  l)d  ■ 


2 kbps 
["8  kbps 


Qd  i 0 
Qd  - 0 
0 < Q < 150 


Qd  > 150. 


Results  on  average  packet  delay  as  a function  of  data  packet 
arrival  rate  for  schemes  v0-v4  are  summarized  in  Table  2.  Table  3 shows 
the  average  bit  rate  assigned  to  a voice  user  for  each  case.  The  first 
scheme  represents  an  improvement  over  no  flow  control  for  packet  arrival 
rates  Ob  ■ 6.0,  7.0,  and  8.0  packets/ frame.  (Recall  that  the  frame 
duration  b ■ 10  ms.)  The  price  for  this  improvement  was  a drop  in 
average  voice  bit  rate  from  8 to  5.7  kbps.  However,  the  flow  control 
worsened  the  situation  when  packet  arrivals  reached  9.0  per  frame.  When 
16  kbps  was  eliminated  as  a voice  rate  (scheme  v2),  a general  improve- 
ment in  delay  performance  was  realized.  In  addition,  the  average  voice 
bit  rate  increased  to  6.6  kbps  because  fewer  users  had  to  be  assigned  to 
the  lower  rates.  Scheme  v3  (the  same  as  vl  except  that  some  data  queue 
monitoring  was  introduced)  represented  a significant  improvement  over  vl 
in  terms  of  data  delay.  Scheme  v4,  where  voice  rate  control  was  based 
strictly  on  data  queue  size,  was  generally  the  most  effective  of  the 
voice  rate  control  techniques. 

However,  none  of  the  performance  results  shown  in  Table  2 are 
really  satisfactory  at  the  high  packet  arrival  rates.  Even  for  scheme 
v4,  the  delay  of  .63  sec  at  9.0  paekuts/frame  is  unacceptable,  and 
corresponds  to  an  average  queue  size  of  571  packets.  Observed  maximum 
queue  sizes  were  significantly  higher.  The  conclusion  is  that,  although 
voice  rate  control  alone  can  enhance  data  performance,  some  form  of  data 
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flow  control  Is  necessary  to  keep  delays  and  queue  sizes  within  reason- 
able limits. 

The  emphasis  here  has  been  on  the  effectiveness  of  voice  rate 
control  in  reducing  data  packet  delay.  However,  the  ability  of  voice 
users  to  operate  at  a variety  of  rates  depending  on  traffic  loads  would 
enhance  the  voice-traf f ic  handling  capability  of  the  system,  independent 
of  its  effect  on  data  traffic.  For  example,  the  assignment  of  lower  bit 
rate  to  new  voice  users  as  the  number  of  active  calls  increases,  could 
be  utilized  as  a means  for  reducing  the  blocking  probability  for  voice 
calls  during  busy  periods. 

B.  Experiments  with  Combined  Data  Flow  Control  and 
Voice  Flow  Control 

The  first  experiment  including  data  flow  control  involved  limiting 
the  data  buffer  to  a fixed  maximum  size  but  not  Including  any 

voice  rate  control.  Qmax  was  (somewhat  arbitrarily)  set  equal  to  150 
packets.  Data  packets  were  denied  entry  to  the  multiplexer  when  its 
queue  was  full.  This  represents  added  delay  since  these  packets  have  to 
be  retransmitted  to  the  multiplexer.  If  the  packets  enter  the  multi- 
plexer directly  from  a user  terminal,  then  the  terminal  would  retransmit 
when  no  acknowledgement  was  received.  Otherwise,  a store-and- forward 
node  feeding  the  multiplexer  could  handle  the  retransmission.  The 
second  experiment  combined  the  same  data  flow  control  procedure  with  the 
data-queue- independent  voice  rate  control  employed  in  scheme  v2  above. 
The  third  test  combined  limitation  of  the  data  buffer  to  150  packets 
with  data-queue-dependent  voice  rate  control  similar  to  the  method  used 
in  scheme  v4  above.  If  the  queue  was  empty,  a rate  of  8 kbps  was 
assigned.  If  the  queue  was  non-empty  but  did  not  exceed  75  packets,  4 


kbps  assigned.  Otherwise,  2 kbps  was  used.  The  three  combined  data/ 
voice  flow  control  techniques  Just  described  will  bo  referred  to  as  dl, 
d2,  and  d3,  and  are  suimnarlzed  as  follows: 

dl,  d2,  d3:  all  Impose  fixed  limit  Q " 150  packets  on  data 
buffer  size.  max 

dl:  R - constant  - 8 kbps 

VCV  £ .5 

•5  < %/Cy  < .75 
.75  < \/Cy  < 1.0 
Qd  - 0 

0 < <!d  i W2 

*W2  < i Qm.x’ 

The  average  packet  delays  for  the  three  cases  are  plotted  in 

Fig.  8.  In  ull  cases,  the  Improvement  In  data  packet  delay  is  dramatic 
as  compared  to  the  situation  where  no  data  flow  control  is  imposed.  For 
example,  the  worst  case  in  Fig.  8 is  an  average  delay  of  49  msec  for 
scheme  dl  at  9.0  packets/frame  arrival  rate.  The  corresponding  best 
result  without  data  flow  control  is  an  average  delay  of  630  msec  for  v4 
as  shown  In  Table  2. 

Tables  4 and  5 show  the  average  voice  bit  rates  and  the  percentages 
of  data  packets  arriving  when  the  queue  is  full,  as  a function  of  packet 
arrival  rate  for  the  three  data  flow  control  schemes.  The  tables 
include  some  measurements  for  packet  arrival  rates  greater  than  9.0 
packets/frame.  Table  5 shows  that  the  percentage  of  packets  finding  a 
full  data  queue  Is  rather  low  in  all  cases.  This  implies  that  the 
additional  delay  seen  by  the  user  who  must  retransmit  these  packets  also 
ought  to  be  modest. 
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new 


d2:  R 


new 


d3:  R 


new 


8 kbps 
4 kbps 
2 kbps 
8 kbps 
4 kbps 
2 kbps 


ii 


DATA  BUFFER  LIMITED  AT  150  PACKETS 


VOICE  RATE 

— — CONSTANT  AT  8kbps  (dl) 

CONTROLLED  BY  VOICE  UTILIZATION  (d2) 

CONTROLLED  BY  DATA  QUEUE  (d3) 


TABLK  4 


If  schemes  dl-d3  are  compared  on  Che  basis  of  data  packet  delay 
alone,  ic  is  clear  from  Fig.  8 that  the  best  performance  is  achieved 


with  d2,  which  includes  voice  rate  control  independent  of  the  data 
queue.  However,  the  improvement  in  delay  performance  over  d3  is  counter- 
balanced by  the  fact  that  higher  average  voice  bit  rates  were  assigned 
in  scheme  d3  at  all  data  packet  arrival  rates  tested. 

Thus  far,  performance  of  the  three  data  flow  control  schemes  has 
been  displayed  on  the  basis  of  packet  arrival  rate  for  convenient  com- 
parison with  the  voice  rate  control  schemes.  In  order  to  take  into 
account  the  effects  of  packet  discard  due  to  data  flow  control,  and  to 
provide  a convenient  means  for  comparing  systems  with  different  capa- 
city, it  is  useful  to  present  the  results  as  a function  of  utilization. 

It  is  also  desirable  to  focus  on  utilization  of  the  variable  capacity 
available  to  data  due  to  fluctuations  in  voice  traffic,  rather  than  on 
utilization  of  fixed  channel  capacity  dedicated  to  data.  A corres- 
pondence between  utilization  of  variable  data  capacity  and  packet 
arrival  rate  is  shown  in  Table  6.  Variable  data  capacity  is  defined  as 
the  difference  between  the  total  capacity  allotted  for  voice  (80  kbps) 
and  the  average  capacity  actually  utilized  by  voice.  The  average  packet 
delays  for  the  three  schemes  as  a function  of  variable  data  capacity 
utilization  are  plotted  in  Fig.  9.  It  will  be  noticed  that  scheme  d3 
performs  more  favorably  than  scheme  d2  in  terms  of  packet  delay  as 
utilization  exceeds  70Z.  Similarly,  the  percentages  of  data  packets 
finding  the  queue  full,  which  are  plotted  in  Fig.  10,  show  that  d3  again 
has  the  best  performance  for  higher  utilization  of  variable  data  capa- 
city. However,  there  is  a price  to  be  paid  by  using  scheme  d3.  As 
utilization  is  Increased  by  greater  packet  arrival  rates,  the  data  queue 
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tends  to  remain  more  than  half  full  during  periods  of  higher  voice 
utilization.  Consequently,  more  calls  are  assigned  at  2 kbps;  and,  lest 
the  average  voice  bit  rates  shown  in  Table  4 be  misleading,  it  should  be 
mentioned  that  for  all  packet  arrival  rates  tested  almost  all  calls  were 
assigned  at  either  8 kbps  or  2 kbps  using  scheme  d3.  This  is  in  con- 
trast to  scheme  d2  where  all  calls  were  assigned  at  8 kbps  or  4 kbps 
with  a ratio  of  slightly  less  than  two  to  one,  regardless  of  packet 
arrival  rate.  If  voice  users  could  tolerate  a greater  overall  incidence 
of  2 kbps  voice  communication  with  a better  chance  of  being  assigned  at 
8 kbps,  it  can  be  argued  that  scheme  d3  is  slightly  better  than  scheme 
d2.  The  percentages  of  calls  assigned  at  8 kbps  are  plotted  in  Fig.  11. 
The  area  of  particular  interest  is  that  corresponding  to  a utilization 
of  variable  data  capacity  just  over  70%  which  is  where  scheme  d3  shows 
improvement  over  scheme  d2  in  both  Figs.  9 and  10.  In  general,  com- 
parable performances  could  be  achieved  with  both  schemes  d2  and  d3. 

Scheme  d3  was  tested  for  a larger  population  of  users  by  increasing 
the  channel  capacity  from  120  kbps  to  600  kbps,  with  a voice  capacity  of 
400  kbps  and  a dedicated  data  capacity  of  200  kbps.  The  offered  voice 
traffic  for  the  larger  population  was  set  at  40  Erlangs.  (With  30 
8 kbps  channels  of  voice  capacity  and  40  Erlangs  of  traffic,  the  blocking 
probability  computed  on  the  basis  of  fixed  voice  rate  operation  at 
8 kbps  is  =*  .0186,  approximately  the  same  as  the  blocking  probability 
for  the  10-channel,  3 Erlang  case  considered  up  to  now.)  The  perform- 
ance of  d3  with  this  larger  population  was  compared  to  that  with  the 
smaller  population  considered  above,  where  80  kbps  was  allocated  for 
voice  and  the  offered  voice  traffic  was  5 Erlangs.  As  shown  in  the  plots 
of  average  data  packet  delay,  percentages  of  packets  finding  the  queue 
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full,  and  the  percentages  of  calls  assigned  at  8 kbps  in  Figs.  12,  13, 
and  14,  respectively,  the  performance  was  Improved  for  the  larger 
population.  This  can  be  attributed  in  part  to  the  more  rapid  fluctua- 
tion in  voice  utilization,  which  causes  the  peak  of  voice  utilization  to 
be  shorter  in  duration. 


V. 


CONCLUSIONS 


An  error  In  a previous  analysis  has  been  identified,  and  modified 
results  have  been  presented  for  the  data  packet  delay  performance  of  an 
integrated  voice/data  multiplexer  structure.  The  new  results  were 
developed  both  through  a simulation  and  analysis,  and  showed  that 
excessive  data  queues  and  delays  could  build  up  under  heavy  loading 
conditions.  A variety  of  flow  control  mechanisms  to  prevent  these  large 
delays  were  investigated.  These  mechanisms  included  voice  rate  control, 
limitation  of  the  data  buffer,  and  combinations  of  voice  and  data  flow 
control.  The  best  results  were  obtained  with  limitation  of  the  data 
buffer  combined  with  data-queue-dependent  voice  rate  control. 


41 


APPENDIX  I 


. Error  In  Previous  Analval.8 

The  analysis  aC  issue  was  developed  in  Equations  (9)  — (13)  of  [1]. 
v D 

and  Qr  were  defined  as  the  number  of  occupied  voice  channels  and 
number  of  data  customers,  respectively,  in  the  system  just  after  the  nth 
gate  opening  (where  gate  openings  are  time  instants  occurring  at  b - 10  ms 
intervals  when  data  and  voice  customers  are  admitted  to  the  system) . An 
expression  was  written  for  the  z-transform  of  the  conditional  probability 
distribution 
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and  this  probability  was  "unconditioned"  in  two  stages  according  to 
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- E Pr[<£  - S-k]  E Pr[Q°  - j]  Pr[Qn+1  - i | Q^«j,Q^-S-kl  (AI.l) 
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However  this  unconditioning  procedure  would  be  valid  only  if  and  Qr 
were  statistically  independent.  One  should  not  assume  such  Independence 
to  hold,  particularly  since  large  tends  to  be  correlated  with  large 


The  correct  unconditioning  relation  is 
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but  it  is  difficult  to  proceed  beyond  (AI.2)  since  the  required  joint 

v D 

probability  distribution  of  Qn  and  is  unknown.  In  proceeding 


A-l 


from  (AI.2),  it  was  assumed  that  Pr[Q^  - k)  - Pr[Qv  - k]  for  all  n. 

This  assumption,  togathar  with  tha  incorrect  unconditioning  procedure, 
has  tha  affect  of  treating  as  a random  procasa  which  is  Independent 
between  adjacent  10  msec  frames.  Actually,  as  discuased  above  and 
illustrated  for  example  in  Fig.  2,  changes  only  with  call  Initiation 
or  termination  and  is  vary  highly  correlated  across  frama-length  time 
intervals. 
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In  order  to  develop  an  analytic  model  for  the  SENET  multiplexer 
structure,  as  defined  In  [2],  we  need  to  assume  exponentially  distri- 
buted rather  than  constant  data  packet  lengths  and  Ignore  the  effect  of 
time  quantization  introduced  by  the  frame  structure.  Under  these 
assumptions  the  system  can  be  modelled  as  a two-dimensional  Markov 
chain.  The  solution  then  is  found  by  the  use  of  generating  functions. 
This  method  is  standard  for  such  problems  and  has  been  used  on  problems 
similar  to  the  one  addressed  here.^’^’^’^’®  In  fact  Chang**  has  analyzed 
the  same  system  as  the  one  addressed  here.  We  include  here  a formu- 
lation of  the  two-dimensional  difference  equations  for  the  system  state 
probabilities  and  a discussion  of  the  general  solution  procedure  for  the 
infinite  buffer  case.  An  explicit  solution  is  presented  for  the  single- 
channel case.  Finally,  results  of  direct  numerical  solution  of  the 
system  difference  equations  for  the  finite  buffer  case  are  presented  and 
compared  with  simulation  results. 

The  difference  equations  for  the  system  state  probabilities  are 
formulated  as  follows.  Let  class  1 be  voice  and  class  2 be  data.  We 
assume  that  the  arrival  process  for  each  is  Poisson  with  parameter  X^, 
i - 1,2;  and  that  the  service  distribution  for  each  class  is  exponential 
with  rate  p^.  If  the  service  and  arrival  distributions  are  independent 
and  is  the  steady  state  number  of  class  1 customers  in  the  system, 
let  Pr{Q^  • i,  Q2  • j}  - Pt  j.  There  are  S + N channels  with  S set 
aside  for  the  voice  calls  (class  1),  with  priority.  In  other  words  if  a 
voice  call  shows  up  and  there  are  less  than  S voice  calls  occupying  the 
channels  the  call  will  seize  a free  channel  if  there  is  one,  otherwise 
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It  will  preempt  a data  customer,  who  might  be  using  one  of  them.  The 

preempted  data  customer  returns  to  the  buffer,  and  is  serviced  on  a 

first-come,  first-served.  Arriving  data  customers  can  use  any  free 

channel  for  the  special  case  where  N ■ 0 (see  [7]).  We  assume  that 

there  are  M buffer  spaces  available  for  the  data  customers.  The  steady 

state  equations  for  P.  are  for  i - 0,1,..., S - 1 

*•*  J 


(X1  ♦ x2  ♦ iUj  ♦ Ju2)Pi>:1  - XlPl_  + (i  + Du^., 


+ X2pi,j-1  + + 1)p2pi,j+l:  J " O*1*---.  N + S i 1 


(All. la) 


(X 

+ 


1 + X2  + iPj  + (N  + S - i)P2)PlfJ  - XlPl.ltJ 


(i  + »"lpi+l,i  + X2Pi,j-l  + (N  + S - i^2Pi, j+1: 


(All. lb) 


N + S-  i<j<M  + N + S-  i-  l 


and  for  J"M  + N + S"1 

(Xi  + ih1  + (N  + S - i)U2)PifJ  - X1Pi_1  §J  ♦ (i  ♦ 1)VW|J 

(All. lb’) 


For  i • S,  we  have 


- 


* 


(All. 2a) 


(All. 2b) 


(X2  + SU1  + JV*2)PS  j - X1ps_1>j  + X2PS,j-l 


+ (j  + 1)vj2ps.i+r  i “ 0,1 N " 1 ; 


(X2  + SW1  + NU2)PS,j  “ XlPS-l,j  + X2PS,j-l  + Ny2PS,j+l: 


N<J<M  + N-1; 


and  for  j ■ M + N 


<s“l  + Nl*2)PS.J  " XlPS-l,J  + X2PS,j-l  ' (All. 2b') 


ln  th'“  e’““0“  P-1,J  5 ps+l.j  5 Pl-1  5 pl.(H*M+l-l  ‘ °'  *nd  “he" 

N - 0 equation  (All. 2a)  is  meaningless. 

The  solution  of  this  system  of  equations  can  be  carried  out  formally 
using  transform  methods  for  the  case  where  M - The  procedure  is  to 
define 


P (z)  - E P.  z 

j-0 


J 


i - 0, 1, . . . ,S 


(All. 3) 


and  use  AI1.1  and  All. 2 to  develop  a set  of  linear  equations  of  the  form 
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A(z)  P(z)  - B(z) 


(All. A) 


where  A(z)  is  an  (S  + 1)  x (S  + l)  matrix  depending  on  the  known  system 


parameters  A^,  A^,  N,  and  S;  P(z)  is  an  (S  + 1)  x 1 vector  with 

components  Pq(z),  P^(z),...,  P^,(z);  and  B(z)  is  an  (S  + 1)  x 1 vector 

which  depends  on  the  unknown  quantities  P.  . for  i ■ 0,1,..., S and 

* » J 

J ■ 0,1,... ,S  and  J » 0,1,... ,N  + S.  Equation  (All. la)  can  be  used  to 
evaluate  recursively  these  P.  . so  that  B(z)  depends  only  on  the  N + S 

J 

unknowns,  P_  for  J ■ 0,1, N + S - 1.  To  determine  these  unknowns, 
v»  j 

N + S equations  are  needed.  One  equation  can  be  found  by  the  condition 

a 

that  the  carried  load  must  equal  the  expected  number  of  busy  servers. 
Another  N - 1 equations  can  be  obtained  from  (All. 2a).  The  remaining 

4 

equations  are  obtained  by  finding  the  S unique  roots  in  (0,1).  Then  we 
require  that  det(A^(z))  also  vanish  at  these  points,  where  A^(z)  is 
obtained  by  replacing  the  ith  column  in  A(z)  with  B(z). 

Once  Pj(z)  has  been  found,  the  expected  number  of  data  customers 
E{Qj}  and  the  expected  waiting  time  E(W^}  are  derived  from 


E{Q2)  - l P (z) 
i-0 


z - 1 


(All. 5) 


e{w21  - \2  e{q2>. 


(All. 6) 


This  solution  procedure,  though  formally  tractable,  is  difficult  to 

U 

carry  out  in  practice  for  reasonably  large  values  of  N and  S.  Chang 
has  worked  out  some  specific  examples  and  developed  several  approxi- 
mation techniques. 

We  will  not  proceed  further  here  with  the  general  infinite  buffer 
problem.  However,  it  is  instructive  to  consider  the  special  case  of 


N • 0 and  S ■ 1.  For  this  case  the  birth  and  death  equations  for  P 
become 


ij 


(^1  + ^2^0,0  " ^1^1,0  + ^2^0,1 


(Xx  + X2  + W2)P0>J  - + U2P0>j+1  + Vo,J-l! 


(All. 7) 


j * 1.2, .. . 


and 


(X2  + U^P^q  - ^^o.o 


(X2  + U^P^j  - Vo.j  + X2Pl,j-l:  j " 1,2,1 


(All. 8) 


Using  these  equations  we  have 


(-X2z  + (X^  + X2  + y2>z  - u2)PQ(z) 


- ^zP^z)  + u2(t  - 1)P0>0 


(All. 9) 


and 


X1PQ(z)  - [X2(l  - z)  + M1]P1(z).  (All. 10) 

The  only  unknown,  Pq  q,  can  be  found  via  the  carried  load  condition  to 

be 
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(All. 11) 


P0,0  “ (1  ‘ p2  " P1P2)/(1  + pl>- 

Returning  to  (All. 9)  and  (All. 10)  to  find  E(Q0)  we  get 

p,(l  + p.)2  + p.X./y 

E(92)  " (1  + pjHl  - pj  - PjP2)  . (All.  12) 

From  this  expression,  and  more  specifically  in  the  factor  X2/|Ji,  one  can 
see  the  dependence  of  E(Q2)  on  the  voice  holding  time  1/y^. 

The  main  analysis  tool  that  was  used  in  this  paper  was  the  simula- 
tion model  discussed  in  Section  II.  In  Section  III,  a numerical  comparison 
with  the  analytic  results  for  the  special  case  S = 1,  N = 0 was  made. 

We  close  this  appendix  with  a comparison  of  the  simulation  results  and  a 
numerical  solution  of  the  birth  and  death  equation  for  the  case  of  S * 

5,  S » 10  and  M *>  150.  This  numerical  solution  was  carried  out  by 
direct  iteration  of  the  two-dimensional  difference  equations  presented 
above  for  the  case  of  finite  M,  together  with  application  of  the  con- 
dition that  the  state  probabilities  must  sum  to  1.  The  results  are 
given  in  Figs.  15  and  16,  and  one  can  see  there  is  reasonably  good 
agreement  between  the  simulation  and  numerical  results.  Exact  agreement 
is  not  expected  because  of  differences  in  the  models  used.  In  parti- 
cular, we  recall  that  exponentially-distributed  packet  lengths  were 
assumed  in  the  analytic  model  and  f ixed  packet  lengths  were  used  in  the 
simulation. 
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Fig. 15.  Comparison  of  simulation  results  with  numerical  solution  of 
birth-death  equations  for  mean  packet  delay  with  a finite  buffer;  S - 10, 

N » 5,  X ■ .05  sec-*-,  data  queue  limited  to  M ■ 150  packets.  The 
numerical  solution  assumes  exponentially-distributed  packet  lengths  of 
average  size  80  bits,  while  the  simulation  used  fixed-length  80-bit  packets. 
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Fig. 16.  Comparison  of  simulation  results  with  numerical  solution  of 
birth-death  equations  for  percentage  of  data  packets  arriving  when  queue 
is  full.  System  parameters  are  the  same  as  in  Fig. 15. 
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